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We consider the effect of the photon radiative correction on the vacuum energy in a superstrong 
magnetic field. The notion of a photon anomalous magnetic moment is analyzed and its connec- 
tion with the quasiparticle character of the electromagnetic radiation is established. In the infrared 
domain the magnetic moment turns out to be a vector with two orthogonal components in correspon- 
dence with the cylindrical symmetry imposed by the external field. The possibility of defining such 
quantity in the high energy limit is studied as well. Its existence suggests that the electromagnetic 
radiation is a source of magnetization to the whole vacuum and thus its electron-positron zero-point 
energy is slightly modified. The corresponding contribution to the vacuum magnetization density 
is determined by considering the individual contribution of each vacuum polarization eigenmode in 
the Euler-Heisenberg Lagrangian. A paramagnetic response is found in one of them, whereas the 
remaining ones are diamagnetic. Additional issues concerning the transverse pressures are analyzed. 
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I. INTRODUCTION 

Large magnetic fields |B| 3> B c , B c = m 2 /c = 4.42 • 
10 13 G (m and e are the electron mass and charge, respec- 
tively) in the surface of stellar objects identified as neu- 
tron stars [U-HI might provide physical scenarios where 
quantum processes predicted in such a regime could be- 
come relevant for astrophysics and cosmology. According 
to quantum electrodynamics(QED) in strong background 
fields, the most important effects are likely to be pair cre- 
ation photon splitting and photon capture [§- 
[l3j . The last two essentially depend on the drastic depar- 
ture of the photon dispersion relation from the light cone 
curve, due to the vacuum polarization tensor n^„ which 
depends on both the Landau levels of virtual electron- 
positron pair, as well as on the external magnetic field. 
As a result, the issue of light propagation in empty space, 
in the presence of B, is similar to the dispersion of light 
in an anisotropic "medium" , with the preferred direction 
corresponding to the external field axis. The phenomeno- 
logical aspects associated with this issue has been studied 
for a long time and recently, the effects of the vacuum po- 
larization tensor on the Coulomb potential in superstrong 
magnetic field have been considered as well |L5 - llq ] . How- 
ever, the problem concerning the magnetism carried by 
a photon has not attracted sufficient attention, except- 
ing [17j . where an photon anomalous magnetic moment 
has been pointed out. The authors of this reference at- 
tempted to derive this quantity in two different regimes 
of the vacuum polarization tensor. On the one hand for 
low energies in weak fields (|B| <c B c ) and on the other 
hand (originally studied in Ref. [HI) near the first pair 
creation threshold and for a moderate fields (|B| ~ B c ). 
However, it was ignored that the photon magnetic mo- 
ment m 1 actually is a vector defined as the coefficient 
of B when the energy is linearly approximated in term 



of the external magnetic field. In consequence, neither 
any discussion about the connection between this quan- 
tity and its angular momentum was presented nor any 
comment about its precession around the external field 
axis was made. 

In the present work we reveal the asymptotic condi- 
tions where the concept of a photon anomalous magnetic 
moment may be adequate. We will analyze the cases of 
superstrong magnetic fields (|B| 3> B c ) for both the low 
and the high energy limit of the vacuum polarization ten- 
sor. In these domains, one eigenvalue of n M „ depends lin- 
early on the external field. In consequence, the Maxwell 
equations in the "medium" seem to describe a massless 
particle with a magnetic moment. In this case m 7 may 
become a characteristic quantity since the external field 
generates a Lorentz symmetry break down which induces 
the nonconservation of the photon helicity. We shall see 
that this notion equips us with an intuitive tool to un- 
derstand, in a phenomenological way, the quasiparticle 
behavior of a photon propagation mode in a highly mag- 
netized vacuum. 

The second purpose of this work is to show that vir- 
tual photons are a source of magnetization to the whole 
vacuum. We will address the question in which way the 
virtual electromagnetic radiation contributes to a mea- 
surement of the vacuum magnetization and therefore to 
increase the external field strength. This aspect might 
be important for astrophysics since the origin and evolu- 
tion of magnetic fields in compact stellar objects remains 
poorly understood [l9j . Some investigations in this area 
provide theoretical evidences that |B| is self-consistent 
due to the Bose-Einstein condensation of charged and 
neutral boson gases in a superstrong magnetic field [20l — 
[23| . In this context, the nonlinear QED-vacuum pos- 
sesses the properties of a paramagnetic medium and 
seems to play an important role within the process of 
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magnetization in the stars. Its properties have been stud- 
ied also in [26T - |29| for weak (|B| <C B c ) and moderate 
fields (|B| ~ B c ) in one-loop approximation of the Euler- 
Heisenberg Lagrangian New corrections emerge by 
considering the two-loop term of this effective Lagrangian 
[30H341 ] which contains the contribution of virtual pho- 
tons created and annihilated spontaneously in the vac- 
uum and interacting with the external field through the 
vacuum polarization tensor. 



The two-loop term of the Euler-Heisenber g L agragian 
was computed many years ago by Ritus [30, 31]. A few 
years latter, Dittrich and Reuter [32j obtained a sim- 
pler integral representation of this term and showed that 
their results agreed with those determined by Ritus in the 
strong magnetic field limit. In the last few years, it has 
been recalculated by several authors using the worldlinc 
formalism [H, [13] and it has been extended to the case of 
finite temperature as well [35j]. Nevertheless, in all these 
works it is really cumbersome to discern the individual 
contributions given by each virtual photon propagation 
mode to the Euler-Heisenberg Lagrangian which should 
allow to determine the global magnetic character gener- 
ated by each of them. In this paper we compute these 
contributions separately. In addition we derive the indi- 
vidual vacuum magnetizations provided by each virtual 
photon propagation mode for very large magnetic fields 
(|B| 3> B c ). Whilst the complete two-loop contribution 
is purely paramagnetic we find a diamagnetic response 
coming from that n^„— eigenvalue which is used to ob- 
tain the photon anomalous magnetic moment. 



In order to expose our results we have structured our 
paper as follows. In Sect. II we recall some basic fea- 
tures of photon propagation in an external magnetic field. 
Some aspects of the Noether currents associated with the 
spacetime symmetries in presence of an B are discussed in 
Sect. III. The corresponding calculations are carried out 
in a basis that diagonalizes the polarization tensor. The 
posibility to define a photon anomalous magnetic mo- 
ment in an asymptotically large magnetic field is studied 
in Sect. IV. We shall show that m 7 turns out to be a vec- 
tor with two orthogonal components in correspondence 
with the cylindrical symmetry imposed by the external 
field. Moreover, we will see that a nonvanishing torque 
exerted by B might generates a precession of m 7 which is 
a manifestation of the reduction of the rotation symmetry 
and, consequently Lorentz symmetry. In Sect. V we per- 
form the calculation of the two-loop contributions to the 
Euler-Heisenberg Lagrangian given by each virtual mode. 
From them we obtain the modified vacuum energy for 
very large magnetic fields in Sect. VI. The correspond- 
ing vacuum magnetizations and magnetic susceptibilities 
are analyzed as well as the transversal pressure. There 
is additional discussion given in the conclusions while es- 
sential steps of many calculations have been deferred to 
the appendices. 



II. GENERAL REMARKS 

A. Symmetry reduction and diagonal 
decomposition of vacuum polarization tensor 

The standard relativistic description of a photon re- 
flects the underlying symmetry of the Poincare group 
(ISO(3, 1)). Its irreducible representations are character- 
ized by two observable degrees of freedom, corresponding 
to the helicity values A = J • k/|k| = ±1. Here, J encodes 
the generators of the 3— dimensional rotation group and k 
is the photon momentum. In the presence of an external 
electromagnetic field, space-time is no longer isotropic 
and Poincare invariance breaks down. In a reference 
frame where the field is purely magnetic, the symmetry 
breaking for neutral particles has the following pattern: 



ISO(3, 1) + B -> ISO(2) x ISO(l, 1). 



(1) 



The effective physical configuration of Minkowski 
space manifests itself as the direct product of the 
2— dimensional Euclidean group ISO(2) and the (1 + 
1)— dimensional pseudoeuclidean group ISO(l, 1). The 
two resulting symmetry groups are associated with the 
transversal and pseudoparallel planes with respect to 
B— direction. The conserved quantities of any uncharged 
relativistic quantum field within this reference frame are 
related to the Casimir invariants of this direct product of 
groups [361 ]. 

For the electromagnetic radiation this type symme- 
try reduction occurs as soon as radiative corrections 
are considered. In such a case an observable pho- 
ton interacts with the external magnetic field through 
the virtual electron-positron (e^) pair whose Green's 
functions S^(x, x'|B) determine the vacuum polarization 
tensor n M „(x,x'|B) = -c 2 Tr [ 7m S?(x, x'|B) 7l/ Sf(x', x|B)] . 
Consequently, the photon behaves like a quasiparticle 
embodying both radiation and properties, the latter 
quantized by the Landau levels. 

The diagonalization of is expressed as 
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with its renormalized eigenvalues Xi depending on the 
scalars, z x = kF* 2 k/(2A) and z 2 = -kF 2 k/(2A), which 
together with k 2 = zi + z 2 form the Casimir invariants 
of the ISO(2) x ISO(l,l) Lie algebra. Here F» v * = 
l/2& vpa ¥ p<J represents the dual of F^,,, whereas A = 
l/4F Ml/ F Miy = 1/2|B| 2 is one of the external field invari- 
ants (the remaining one, $ = l/4F /iI/ F* t/ , vanishes iden- 
tically). 

In one-loop approximation, the eigenvalues of read 
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Here and in the following s = eBr, 

1 sinh(s) cosh(s7/) — ry sinh(s?/) cosh(s 



(4) 



o-i (s, 77) 



CT 2 (S, ?7) = 



sinh(s) 



■cosh(s), 0-3(8,77) = — 



M(s,77) 



. cosh(s) — cosh(s??) 

Ms, 77 = V- u \ > N M 

2 smh(s) 



sinh(s) ' 
1 - rf 



(5) 



The four vector aft in Eg . ([2]) denotes the correspond- 



ing eigenvectors of [37j : 
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k 2 F 2 A k* - (kF 2 k)k M 
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which fulfill both the orthogonality condition: 

5 1 * (a^) and the completeness relation: 8^ v = 

Eti^ (i) a«/(a«) 2 . 

Owing to the transversality property (k^n^ = 0), the 
eigenvalue corresponding to the fourth eigenvector van- 
ishes identically (x^ = 0). Consequently, the photon 
propagator can be decomposed as 



(XV 



3 1 



a (i) a (i) 



k 2 -x« ( a (0)< 



k 2 k 2 



(7) 



with £ being the gauge parameter. 

According to Eq. ([7]) three nontrivial dispersion rela- 
tions arise 



.XT, (z 2 , zi, A) for i = 1,2, 3. 



(8) 



We should keep in mind that the general structures of 
Z2 and zi are complicated (see Appendix |X| since both 
of them depend on the relative B— orientation with re- 
spect to a reference frame. Substantial simplifications 
are achieved in reference frames which are either at rest 
or moving parallel to the external field. In these cases 
Z2 = k 2 ± and zi = k 2 — lu 2 . Here k^ and k|| are the com- 
ponents of k perpendicular and along the external field 
respectively, with k 2 = k^ +k 2 — lu 2 . Note that, as a con- 
sequence, x, depend on both A = 1/2(B 2 + B 2 +B 2 ) and 
on the B-direction with respect to our reference frame. 
Solving Eq. © for k 2 = k 2 — uj 2 in terms of Z2 yields 



k 2 +m 2 (z 2 ,A), 



(9) 



with the term arising as a sort of dynamical mass. 
Note that vanishes for z 2 = due to the gauge invari- 
ance condition ^(0,0, A) = 0. Obviously, the dispersion 
law given in Eq. © differs from the usual light cone 
equation. This difference increases near the free pair cre- 
ation thresholds remarking the quasiparticle feature of a 
photon in an external magnetic field. For more details 
we refer the reader to Ref. 

By considering (k) as the electromagnetic four vec- 
tor describing the eigenmodes, we obtain the correspond- 
ing electric and magnetic fields of each mode 

e«=j(u (l 'a< i »-kaf) and h^=-ikxa«. (10) 

Up to a non essential proportionality factor, they are 
explicitly given by: 
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= — in±u>, = iky x n^, 

iki_k||, e^ 2 ' =m||(k| -cj 2 ), 

h/ 2 * 1 = iu(k± x nil), 

e^ 3 ) = icj(n^ x nil), 



(11) 



ini|k_L, h 
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-inj_k||. 



Here, ri|| = k||/|k||| and = kj_/|k_i_| are the unit vec- 
tors associated with the parallel and perpendicular direc- 
tion with respect B. 



III. NOETHER CURRENTS OF THE 
RADIATION FIELD IN PRESENCE OF B 

A. The Poynting vector and the physical 
propagation modes 

Whenever the Minkowski space is translations in- 
variant, the associated Noether current of the electro- 
magnetic radiation (A M (x)) provides a conserved stress- 
energy tensor: 
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TV 



(12) 



Here = 9 M A„ — c^A^ is the electromagnetic tensor, 
whereas the metric tensor tf ,v has signature + + H — with 
77 11 = rj 22 = 7/ 33 = — rf° = 1. In this context the momen- 
tum density is defined by ^ = T° M . Note that the sym- 
metry reduction by the external magnetic field does not 
alter the translational group involved within ISO(3,l). 
Therefore, all — components are conserved. 

The spatial part of *}} is the Poynting vector density 



V = E x H. 



(13) 



For each eigenmode it reads J>W = eW x In partic- 
ular 



p«=wk,|, p< 3 )=wk, 



jf^kik,,, p<# = W (o; 2 -kf)kx. 



(14) 



4 



Different photon degrees of freedom contribute in the 
presence of B, depending on the direction of wave prop- 
agation: for a pure longitudinal propagation k || B, the 
Poynting vector p^ = 0. As a consequence, eigenmode 
2 is a pure longitudinal non physical electric wave and 
does not carry energy. On the other hand, the first 
and third mode have well-defined Poynting vectors along 
the external field p^ = p( 3 \ In this case, each set 

jeWjhW^Wj and |e( 3 ) , h (3) , p^ j forms an orthog- 
onal set of vectors and represent waves polarized in the 
transverse plane to B. Consequently, for pure parallel 



propagation and represent physical waves. 

Now, if the photon propagation involves a nonvanish- 
ing transversal momentum component k^ =^ 0, we are 
allowed to perform the analysis in a Lorentz frame that 
does not change the value k_L, but gives k|| = and 
does not introduce an external electric field. As a conse- 
quence, the energy flux of the first eigenmode p^ = 
and becomes purely electric longitudinal and a non phys- 
ical mode. In the same context each set |e' 2 \ h*- 2 -*,^ 2 ) j 

and |e( 3 \ h*- 3 ^, p^ j forms an orthogonal set of vectors 

and represent waves polarized in the transverse plane to 
B. Hence, for a photon whose three- momentum is di- 
rected at any nonzero angle with the external magnetic 

(2) (3) 

field, the two orthogonal polarization states a^ and a^ 
propagate. Note that mode 3 represents a physical wave 
independent of the direction of propagation. 



,(3) 



B. The spin and boost of the eigenwaves 

The Noether current associated with rotational invari- 
ance is usually split into two pieces: orbital angular mo- 
mentum density and intrinsic spin density. For an elec- 
tromagnetic field, the latter can be derived from the third 
rank tensor 



{f a ?)lK{x) (15) 



with Jz? = being the free part of the QED- 

Lagrangian. Here {J a $Y v = -i fe?7/3i/ - ^v) de- 
notes the four-dimensional representation of the Lorentz 
generators. By considering the previous definition we can 
express Eq. (IT5l) as 



the classical spin density of the electromagnetic field 

5(x,t) = -iE(x,t)xA(x,t). (18) 

For /i = (3 = we obtain the corresponding Noether 
current related to the boost transformations, which is 
given by &(x, t) = iE(x, t)A (x, t). 

In order to analyze the behavior of <S and & in pres- 
ence of B we regard (see Eq. (fT0|) ) as the electric 
field associated with the fourth potential (aW ) . Up to a 
nonessential factor of proportionality we have 

*« = -ie« x = -iatfh® and = ia^e^(.19) 



Manifestly Eq. (fT9|) expresses the connection between 
the Noether currents and the different polarization planes 
associated with each eigemode. In particular 



M - 



gr(k|| x k x ), - ^k ± , 



.(2) 



5^ = oj(k\\ x kj 



k i k„ , k 



(2) 



Wk - w 2 



(20) 



5 ( 3 ) = 0, k( 3 ) = 0. 

According to the reduced Lorentz symmetry, i.e. Eq. 
([T]), only the parallel components sjp and k$ are related 
to conserved quantities and therefore, just the component 
of the electric and magnetic field along B can generate 
conserved charges. However, for a purely parallel propa- 
gation (k± = 0), this connection is rather obscure due to 

the absence of e [| 1 an d f" | 1 ' 3 ^ - I n this case the first and 
third mode may be combined to form a circularly polar- 
ized transversal wave which is allowed by the degeneracy 
property: 



xri (z^O,^) = x 3 (zi,0,$). 



(21) 



In this context the photon is labeled by the helicity 
A = J|| - n 1 1 which seems to be a conserved quantity. For 
nonvanishing transversal propagation kx ^ 0, however, 
A stops being a well-defined quantum number due to the 
nonconserved rotations transversal to B. 



C. Connection between helicity and spin for pure 
parallel propagation 



(16) 



We fix /j = in Eq. p^|) and take into account the spa- 
tial part of the remaining tensor. Under such a condition 
the intrinsic angular momentum density of the electro- 
magnetic field is reduced to 



5° =-i \E i (x)A j (x)-E j {x)A i (x)], 



(17) 



with E l (a;) = y° being the electric field. The spatial 
components defined from |e yfc iSy with e 123 = 1 define 



In order to establish the connection between helicity 
and classical spin of the electromagnetic field we express 
S£p =S§ + Here 



{ffg ± iTc?} Kp-\ {r p ± iTf} A a (22) 



and is the dual of Adopting a similar procedure 
to those developed below Eq. (fT6]l the vectors read 



(23) 



S ± = — - {E ± iH} x A. 
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The complex fields | {E ± iH} transform irreducibly un- 
der spin (1,0) and (0,1) representation of SO(3, 1) ~ 
SU(2) x SU(2), respectively. These field combinations 
fulfill the free Maxwell equation for left- and right circu- 
larly polarized radiation 

V x (E±iH) =F^(E±iH) = (24) 

corresponding to A = =pl. 

For propagation purely parallel to the external mag- 
netic field we define the electric field e( c ) = eW + e' 3 ) 
and magnetic field = — ik x a( c ) = tr 1 -' + asso- 
ciated with a' c ) = + sS 3 \ respectively. A this point 
it is meaningful to analyze 

s (c) = -te< G > x a<°> and = x h (c) . (25) 

Inserting the explicit expression of e( c ) , h*- * 1 and aS^ in 
Eq. (H5l) we get 



,W=,(i) + ^(n || k a .-n J .k || ) ) 
+p(3) -u(k ± x n,i). 



(26) 



Note that 



lim s (c) = and lim p (c) = 2uk». (27) 

k_i_->-0 kx->0 " 

The resulting limits are expected for a circularly polar- 
ized wave. 

Further considering Eq. (|23[) we obtain 



D. Spin density of a free electromagnetic field for 
perpendicular propagation 

Let us consider a free photon propagating perpendic- 
ular to the external magnetic field (w(2,3) = |k|). In this 
case, the behavior of the photon spin density reads 

s (t) = _ ie (t) x a (t) = _ ia « h (t) (30 ) 

Here eM = e< 2 ) + e< 3 ) and h (t) = h (2) + h (3) denotes the 
electric and magnetic field associated with aft = a^ + 

a^ 3 \ respectively. Because of a^ = — ky, = and 
Eq. (HI]), we find 



s (t) =s (2)_ nxs (2) 



(31) 

with n = k/|k| being the corresponding wave vector. 

Additionally, taking Eq. (|27|) into account we intro- 
duce the vector 







for k i =0 



s (2) - n x s (2) for k ± ^ 



(32) 



with s^ 2 ) =s^ 2 V|s^ 2 ^| = (nil xrijj. The above expression 
is total spin density, which itself depends on the direction 
of propagation. 

Because of the fact that s^ 2 ) • k = 0, the total spin den- 
sity of the electromagnetic field is orthogonal to the wave 
vector (s 7 -n = 0). We remark that both s' 2 ) and n x 
are orthogonal to each other. Note that for transversal 
propagation the total spin density of the electromagnetic 
field has a parallel component to B given by s| = — j^ n \\ ■ 



IV. PHOTON MAGNETIC MOMENT 



-i(eM±ih< c >) xaW 
-s^ ± -h< c > x a( c > 



Note that 



k'» 



(28) 



h (c »xa< c ' = ik + i kn+i-W 

kr " kr 

- ( k ^ Xn ll)' 

The above relation was obtained by inserting the explicit 
form of e( c \ and a( c ). Therefore, the angular momen- 
tum associated with left and right circular polarization 
are represented by 



lim s± ~ Tkii 



(29) 



In this section we explore the possibility that a pho- 
ton might carry a magnetic moment. We will restrict 
ourselves to asymptotically large magnetic field b = 
B|/B c 3> 1. In this limit the eigenvalue of the second 
propagation mode contains a term linearly growing with 
the magnetic field strength at low and high energy lim- 
its. Both cases deserve a separate study because the 
>t2~- structures differs from one to the other energy do- 
main. Under the same conditions k\ and k% cannot cre- 
ate a virtual electron-positron pair in the ground state 
[TT| . These show logarithmic dependences on |B| and 
their corresponding physical dispersion laws are inde- 
pendent of the external field (for details see Ref. [l3[ ). 
Therefore the cases concerning to the first and third prop- 
agation mode are not relevant in the current context. So, 
in this section, we will analyze the effects produced due 
to the second eigenmode. 



Regarding the normalized version of the above limit 
(s± = T n ||) a s the angular momentum used to define 
the helicity, we can write A = s± • p( c ) = =pl, where 
p( c ' = nil is the normalized version of the second limit 
computed in Eq. (|2"71) . 



A. Infrared struture of xi: Covariant 
decomposition of the photon interaction energy 

In the limit m 2 b ;§> m 2 ;§> u; 2 — k 2 with m 2 b >• k\, 
the second eigenvalues of n MJ/ shows a linear function on 
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the external field strength 



4* 



-£>(b)zi 



« F - k "<''" 2, ■ < 33 » 



with g(b) = -^b and a = e 2 /An = 1/137 being the fine- 
structure constant. Note that we have used the decom- 
position: zi = k"F* A F£"k„/(2A) = k^F* A af '"/(2A) 1 / 2 . 
The expression for is analogous to that of the invari- 
ant interaction energy of the electron [2^ |: 



Here g = ^ is a kind of Lande factor, whereas /[kj_] = 

kj_/m. In addition, S Mt/ = S^/dkl) is tensor with 

fiffi = — ik^al?' + ik^alf-* referring to the antisymmetric 
electromagnetic tensor generated by the third propaga- 

(3) 

tion mode. Note that = e 
The spatial part of 



'i^v can be written as: 



u(3) 

2m k 



(39) 



2m : 



(34) 



Here is the electron fourth momentum and S M = 
7 5 (7^ — — ) is the electron spin with n M = p^ + ^eF^x". 
In the electron rest frame, e describes the interaction en- 
ergy between the electron magnetic moment and the ex- 
ternal magnetic field. In the photon case the structure 

of Eq. (|33|) shows that the electron spin S M is replaced 

(2) 

by the photon polarization , describing the intrinsic 
rotation of the photon in the plane perpendicular to the 
external field. Note, however, that the absence of a pho- 
ton rest frame prevents the definition of a photon spin, 
unlike the electron case. To avoid this problem and for 
further convenience we will analyze the photon propaga- 
tion by investigating its dispersion law Eq. ([8]), which 
reads 



g(b) (1 + 0(b))- 1 z 3 . 



(35) 



For magnetic field strength 10 < b <C 3ir/a one should 
treat 0(b) as small. The expansion of w up to first order 
in 0(b) gives the dispersion law 
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U.--2 



0(b) Z2 /|k| 



(36) 



Obviously, the first term in Eq. (|36|) corresponds to 
the light cone equation whereas the second arises due to 
the dipole moment contribution of the virtual electron- 
positron pair. In this approximation the dispersion law 
does not essentially deviate from its vacuum shape and 
its interacting term grows linearly with the external mag- 
netic field. This fact attracts our attention because it 
seems that a magnetic moment may be ascribed to a 
mode-2 photon. 

Hereafter we will denote the interaction energy by 



W = - 



1 0(b)z 2 



a e 



(z 2 ) 



1/2 



67r m 2 



F^k^ 3 > 



WF A k A /(2A) 



where the decomposition z 2 
(z 2 ) 1/2 k^F A1I ,a I/ ( 3 V(2A) 1 / 2 has been used. Note that 
(z 2 ) 1 / 2 = kj^ in reference frames which are at rest or 
moving parallel to B. With this in mind, we find by sym- 
metrization of Eq. ((37j) 
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Jtuv?^ With 



-fc^/ptJS/",. ( 38 ) 



Here ^wl = —ieum 1S the 3— dimensional representa- 
tion of the SO(3) -generators, fulfilling = 
ie iikjtk) and U W )2 = = 2S lm . 

In a Lorentz frame where F Mi , is purely magnetic the 



structure of Eq. 



is expanded to 



7/ 



0^— /[kx] [ri|| sin — n_LCOS(£] • B, (40) 



where < 4> < tt is the polar angle between the wave 
vector n and the external field (tan0 = kj_/k||). The ex- 
pression inside the brackets is a unit vector orthogonal 
to the direction of propagation. We then write the inter- 
action energy as 



?/ = -m 7 • B 



with 



,— /[kj(nxs 



J 2 ) 



(41) 



(42) 



The struture of the interaction energy is similar to a po- 
tential energy of a magnetic dipole in an external mag- 
netic field. Thus, in first approximation (with respect 
0(b)), the second propagation mode in a superstrong 
magnetic field seems to behave as a quasiparticle having 
a magnetic moment m 7 . This behavior occurs in some 
scenarios of condensate matter and has allowed to in- 
troduce the concept of polariton. The latter results from 
strong coupling of electromagnetic waves with a magnetic 
dipole-carrying excitation. 

According to Eq. (|40|) . this magnetic moment is the 
sum of two orthogonal components, as dictated by the 
cylindrical symmetry of our problem: the first one along 
the external field direction 



2ni /[k±]n||Sin< 



(37) whereas the second one is perpendicular to B 



y 2rn L J 



jnj_ cos< 



(43) 



(44) 



Furthermore, both components of m 7 show opposite 
magnetic behavior: while the parallel one is essentially 
paramagnetic (m\ > 0), the perpendicular one is purely 
diamagnetic (rre 7 < 0). They become nonmagnetic for 
propagation along B, (k^ = 0), in whose case the radi- 
ation is, -like in the free Maxwell theory-, insensitive to 



7 



the magnetic field. Note that for kj^ — )• the photon 
anomalous magnetic moment vanishes. This is a direct 
consequence of the gauge invariance of the eigenvalues 
of 11^ . Clearly, the effective magnetic interaction is re- 
lated to mJj which is invariant under rotation around the 
external field. 
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FIG. 1: Photon anomalous magnetic moment drawn against 
the transverse momentum component. Here jiB = e/(2m) = 
9.274 • 10~ 21 erg/G is the Bohr magneton. 



The dependence of the photon magnetic moment on k j_ 
is displayed in Fig. [T]for purely perpendicular propaga- 
tion with respect to B. Note that in this case m 7 • n = 0. 
Within the range of frequencies for which m 7 is defined, 
it is 4 orders of magnitude smaller than the electron 
anomalous magnetic moment (pf — ^r^m) @ but it is 
12 orders of magnitude larger than the neutrino mag- 
netic moment (fi Ue — 1CP 19 ^b) [HI- However, due to 
the astonishing experimental precision with which the 
anomalous magnetic moment of both the electron and 
muon are measured, there is some hope for an experimen- 
tal -probably astronomical- measurement of the photon 
anomalous magnetic moment. 

Note that although the interaction energy Eq. f4"Tj) 
seems to be linearized on the external field and m 7 is 
independent of |B| it turns out to be a vector depend- 
ing on the external field direction. Therefore, if the 
photon anomalous magnetic moment is interpreted as 
/x 7 = —duj/dB, some differences with respect to m 7 are 
expected. However, in Appendix [Al we find that 

H 1 ~ m 7 . 

We would like to bring the attention that this procee- 
dure identifies the magnetic moment only if the disper- 
sion law depends linearly on the external magnetic field 
(for the electron case see Ref. |39l44l| ). Otherwise /x 7 
must be understood as a sort of "photon magnetization" 
rather than a magnetic moment. However the latter no- 
tion equips us with an intuitive tool for qualitative anal- 
ysis of the magnetization generated by a background of 
observable photons. 



B. Discussion 

We have seen that the external magnetic field leaves 
only the rotational symmetry around of B invariant. 
Consequently, rrv\ is an invariant under the rotation 
around the axis of B while m 7 does not. Note, in ad- 
dition, that m 7 does not contribute to the interaction 
energy. However m 7 turns out to be relevant in another 
case: the external field does exert a torque on the mag- 
netic dipole which tends to line up m 7 with B. Indeed, 

r 7 = m 7 xB = xB = 0^-/[kJ|B| cos0s (2) . (45) 

Because of the fact that the torque is collinear to s^ 2 \ 
the transversal components of total angular momentum 
in that direction is not conserved. This is a manifesta- 
tion of the reduction of the rotation symmetry and, con- 
sequently Lorentz symmetry. Therefore m 7 is a signal 
of a break-down of the Lorentz symmetry. This fact cor- 
roborates and extends the result presented in Ref. [42[ 
which claims that a magnetic dipole moment of truly- 
elementary massive neutral particles is a signal of Lorentz 
symmetry violation. Note, in addition, that r 7 is orthog- 
onal to n x s*- 2 ) which might causes a precession of m 7 
around B in analogy to the Larmor precession for elec- 
trons. 

The Lorentz symmetry breaking implies that the pho- 
ton helicity is no longer a conserved quantity due to ro- 
tations transversal to B (see Subsect. IHI B[) . As a con- 
sequence, alternative conserved elements are needed to 
identify a photon interacting with B. Once the magnetic 
moment is a coefficient of the linear contribution in the 
magnetic field it may become a characteristic quantity 
since one may attribute it an "spin." However this is 
only a possible point of view and only further studies 
can tell how far the interpretation can be stretched. A 
realistic treatment of this problem requires a group theo- 
retical analysis including m 7 similar to those developed 
by Wigner for a free photon (38|. The latter is beyond 
the scope of this manuscript, but a study in this direction 
is in progress. 

Let us consider now a system consisting of ,jV inde- 
pendent photons described by the dispersion equation 

oj = |k| - ra 7 • B (46) 

where we have considered Eq. Eq. ([57)1 and Eq. 

(|41"j) . In addition we will suppose that they propagate 
purely transversal to the external magnetic field. Accord- 
ing to these assumptions, the total interaction energy of 
the system is 

J/ 

«toW = ~X>f } |B| = -^(m«)|B| (47) 

i—i 

(i) 

where m\ denotes the magnetic moment of each pho- 
ton. Here the average of the photon magnetic moment 
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is (m'A 



= E 



I JV ~ (kj_). For a uniform pho- 



ton distribution in the transversal plane to B, one ex- 
pects that (kj_) = which implies that %otai = 0- As 
a consequence, the system does not carry a magnetiza- 
tion 9JT 7 = — <9%otai/<9B. If the system is considered as 
a monochromatic beam (kj_) ^ 0, ^totai 7^ and 



I sot. 



-W to tai/a|B| = .jV{m\) 



(48) 



(/[kj) = (k x >/m 



In consequence the beam carries a nonzero magnetization 
which alters the external field B: 



H = B + 47r(2JI 7 



Here A<f vac > is the vacuum magnetization (for details 
see Sect. IVI[) . Within the magnetic field interval for 
which m 7 is defined, |M vac | - 10 10 - 10 12 erg/(cm 3 G). 
Now, in the frequency range of x-rays ((kj_) > 150 eV) 
the averaged magnetic moment (m 7 ) ~ 10 _27 erg/G. In 
order to produce a photon magnetization of the order 
of ~ |M vac |, a photon density of order p 1 = jY jN ~ 
10 36 - 10 38 

10 42, 



,-3 



cm " would be necessary. Moreover, for 
Pi ^ 10 42 cm~ 3 the magnetization carried by the beam 
is |SDT 7 | ~ 10 15 erg/(cm 3 G) and thus larger than \M V&C \ 
even more, it is 1 order of magnitude larger than the 
external magnetic field |B| ~ 10 14 G/cm 3 . 

Finally, we want to conclude this subsection pointing 
out that Eq. P5|) allows to determine the the corre- 
sponding vacuum refraction index in terms of the photon 
magnetic moment: 



|k| 



rrij ■ B 



(49) 



The latter expression depends on the direction of the 
photon momentum and reaches its maximum for a purely 
transversal propagation (kii = 0) in which case n max ~ 
1 + ^25?r|B|. For b ~ 100 > "max ~ 1-038 which exceeds 
the values of typical gases at atmospheric pressure in 
absence of B [13j. We want to remark that Eq. (|49]) 
might play an important role in the evaluation of the 
magnetic effect on gravitational lenses in the vicinity of 
highly magnetized stellar objects [43|, |44[ . 



Ultraviolet domain 



Next we consider the high-energy regime. In the limit 



m 2 b » ur 



'II 



oo and m 2 b 3> kj_, the second eigen- 



value approaches 



TJV 



2« 

— eB. 

7T 



(50) 



Here tri2 is the photon effective mass [45[ corresponding 
to the topological one in the 2— dimensional Schwinger 
model [46]. Its derivation is closely related to the chi- 
ral limit in which the axial current is not conserved (for 



more details we refer the reader to Ref. [47]). As a con- 
sequence, an ultraviolet photon seems to behave like a 
neutral massive vector boson whose movement is quasi- 
confined in 1 + 1 dimensions 



m: 



1/2 



(51) 



In contrast to the previous case does not depend ex- 
plicitly on the polarization mode. This fact is also man- 
ifest in the dispersion law which cannot be linearized in 
the external magnetic field and therefore a photon mag- 
netic moment is not expected. However, this kind of 
photon carries a magnetization density given by 



Mo 



dB 



a e 



(52) 



Note that /x 7 < behaves diamagnetically and de- 
pends on the external field strength. In particular for 
k 2 <ml 



a e 

AS n || 

7 7T m 2 " 



2a 



1/2 



(53) 



which tends to vanish for |B| — > oo. 

For b ~ 10 3 corresponding to magnetic field |B| ~ 
10 16 G the photon magnetization of an ultraviolet radi- 
ation reaches values of the order of /i 7 ~ — /xb- For the 
same magnetic field strength, an infrared photon prop- 
agating perpendicular to B has fj}\ ~ 0^£>(b)~ 5 / 2 ~ 
10~ 12 zxb (see Appendix IA"| . As a consequence, the mag- 
netic response of a photon background compoused by an 
equal number of infrared an ultraviolet radiations will 
be dominated by the ultraviolet contribution, which for 
mode -2 is diamagnetic. The latter behavior is also man- 
ifest in the contribution of this eigenmode to the vacuum 
magnetization density (see Sect. IVII) . 



V. TWO LOOP TERM OF THE 
EULER-HEISENBERG LAGRANGIAN 

Virtual photons can interact with the external mag- 
netic field by means of the vacuum polarization tensor. 
As a consequence they might be a source of magnetiza- 
tion to the whole vacuum. In what follows we compute 
this contribution for very large magnetic fields (b 3> 1), 
which might exist in stellar objects like neutron stars. 



A. The unrenornamalized contribution due to the 
vacuum polarization eigenvalues 

We start our analysis with the Euler-Heisenberg La- 
grangian 



O _ n(0) , Ml) 



(54) 
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with 



4 0) 



i B 2 

2 



and B = BnZ. 



-1/2 



0^3(lZoop)- 



(55) 



Here is the free renormalized Maxwell Lagrangian in 
a Lorentz frame where the electric field vanishes, E = 0. 
Hereafter B , m and eo will be referred as the "bare 
magnetic field strength" , "bare electron mass" and "bare 
charge" , respectively. Without the index these quanti- 
ties must be understood as renormalized. On the other 
hand, is the regularized one-loop contribution due 
to the virtual electron-positron pairs created and annihi- 
lated spontaneously in vacuum and interacting with the 
external field [5j: 




EH' 



FIG. 2: Two-loop expansion of the Euler-Heisenberg La- 
grangian. The solid lines represent the electron-positron 
Green's functions, whereas the wavy line refers to the pho- 
ton. Here is represented by the one-loop graph which 
gives the contribution of the virtual free electron-positron 
pairs created and annihilated spontaneously in vacuum and 
interacting with the external field. The radiative corrections 
(involved in & 2 ') emerge from the two-loop graph due to ex- 
change of virtual photons. 



a(l) _ 
R 



1 



dr 



scoth(s) - 1 - — ) (56) 



1/2 

with s = cBr and c = CqZ 3 . In this context, the one- 
loop renormalization constant is given by 



Z„,\, s = lim < 1 H In ( ^ — 

HUoop) To ^ \ 37r \ im %T 



(57) 



where 111(7) = 0.577 ... is the Euler constant. 

The contribution due to virtual photons interacting 
with the external field by means of the vacuum polar- 
ization tensor is expressed as 

£(2) = \ / (04 n ^( k )®^( k )- (58) 

Substitution of Eq. ([2]) and Eq. (0 into the latter ex- 
pression gives 



x(2) 



E 



d 4 k 



2^7 (27T)4(k 2 -Xi)* 



(59) 



Manifestly, this expression shows that each photon prop- 
agation mode contributes independently. The quantity 



n^(k)D^^k) = £i^i(k 2 



1 represents the inter- 



(2) 




FIG. 3: Diagrammatic decomposition of £' 2 ' in terms of the 
vacuum polarization eigenmodes. 



into three diagrams (see Fig. [3]). Since we are inter- 
ested in the magnetic properties generated by each pho- 

(2) 

ton propagation mode, each term £ 2 will be studied 
separately. 

The symmetry of our problem suggests to perform the 
k integration with two sets of polar coordinates corre- 
sponding to the two planes (k x ,k y ) and (kn , k^), with 
k4 = iuj. However, it is convenient to use the integration 
variables z 2 = k^ and zi = k 2 + k| > 0, such that 



0(2) = 



32tt 2 J 



d^dz2dzii<i(z2, zi)e 



(Z2+Z1) 



(61) 



Substitution of Sq in Eq. (I61[) yields: 



i(2) _ "0 



cation energy between the full photon propagator with 
the vacuum polarization tensor. In order to obtain the 
leading term in an expansion in powers of e 2 , i. e. the with s = eoBor and 
two- loop contribution as shown in Fig. ([2]), we just ne- 
glect Ki in the denominator of the photon propagator. Qi(s rf) 
Hereafter we consider this approximation, in which case 



•—^3 / drexp(-m 2 r) / drjQ^r,) (62) 
6llT Jo Jo 



sinh s 



{V(s,7y)+W(s,77)}. 



(63) 



£( 2 >=5>f with £f } ^ 



d 4 k 



2 J (2tt) 4 k 2 ' 



(60) 



Note that Xj denotes the unrenormalized eigenvalue. Its 

expression is given by Eqs. (j3][5]) but considering only 

the first line of Eq. Q where all physical parameters are 

unrenormalized. 

Since the analytical properties of any differs from 

(2) 

the others, each contribution £) ; will be different and 
therefore, the original 2— loop graph can be decomposed 



Q2M = ^V(s,7 ? ) + -^-W(s, ?7 ), (64) 
sinh s sinh s 

GsM = -^-V(s,7 ] ) + -^W(s,? 7 ). (65) 
smn s smn s 

The functions V(s, 77) and W(s, 77) are given by the follow- 
ing integral representations 



V(s,77) 



W(s,77) 



1 

d£dz 2 dz 1 z 1 e- Z2 (^ +£ )- Zl (5+^) ! (66) 



/ dMz 2 dz lZ2 e- Z2 (^+ £ )- Zl (^^). (67) 
Jo 
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These integrals can be explicitly calculated: 

V(S,T?) 



28888 *f£), (68) 



W(s,r?) 



(M - N)N (M - N) 2 V N 

2e 3 B 3 , 2e 3 B 3 ^/M, (()()) 



(N-M)M (M-N) 2 V N . 
such that 

V(s, 77) + W(s, 77) = 2e 3 ;B 3 ) M- 1 (s, 7 7 )N- 1 (s, 77). (70) 



To regularize the remaining integration over 7/ we ex- 
press Eq. (|73[) in the following way 



,(2) 



an 2 P° 

-\ / dre- m o- / d7,Q i2 (s,77) (74) 

6Z7r J2t Jo 



32? 
32tt 3 



2t„ 



dre" m o T / d7 ? Joi(s, 77) 



6A7T Jo Jo 



B. Renormalized contributions due to the photon 
polarization modes 

While the r— integral in Eq. (1621) does not diverge at 

t — > 00, the integrand is singular for r — > 0. In order to 

(2) 

regularize £J we introduce a finite lower limit tq > 
for the proper time integral of the electron propagators 
contained in the polarization tensor. As a consequence, 
we can write 



(2) _ gO 



£ i; = - 



32tt 3 



/2t JO 

with 770 = 1 — 2r /r and 



dre- m « r / dT? { Qi (s, 77) - Q i0 (r, 77)} 



(71) 



Qr0(s,v) 



8(e B ) 3 
(1 — ?7 2 )s 3 



(72) 



(2) 

The substraction of this term guarantees that £) ; = 
for vanishing magnetic fields. Obviously, Eq. (I7T1) differs 
from the original two-loop contributions in a term which 
is magnetic field independent and therefore a constant. 

We proceed by adding and subtracting the functions 
Qi2(s, 77) in the integrand Eq. (fTTj) . such that 

(2) _ a 



where Jo^ (s, 77) is a function determined by the singular 
part of the Laurent series of Qj(s, 77) — Qj(s, 77) at r\ = 1. 
As a consequence, the integrations in the third line of Eq. 
(|74"|) converges for t — ► 0. In particular, 



£i(s,77) = 



%(s,r?) 



#3(s,»?) 



4(e B ) 3 (.,(•,, s 2 

^— — scoth(s) H 5— - 2 , (75) 

s 3 (l - 77 2 ) V sinh 2 (s) ' 



2(e B ) 3 
s 3 (l - 77 2 ) 

2(e B ) 3 
s 3 (l- 77 2 ) 



3s coth(s) 



s coth(s) 



sinh (s) 3 



¥-4 1. 



3s 2 2s 2 



sinh (s) 3 



4 



We consider the integration over 77 in the first and sec- 
ond line of Eq. (|74p to express (after some algebraical 
manipulations) 



07T Z J 2 Tn T 3 



dm 2 , 



4 2) 



.^p_ £ (0) 

6w 2 



'2t 

ln( 7 m 2 T ) f°° —e'^ 

J2t t 



r°° dr 2 

J2r T 



Oo 

2?r 



In (7moT ) £^ 



(1) 



-^3 / dre-or / d?/ Q i2 (s,77) (73) + 6° 5m 2 

J2r Jo o 



(77) 



327 
32tt 3 



2r„ 



dre- m 5 r / d7 ? {Q i (s,77)-Q i (s,7?)} 



£l 2) = - 



12tt 2 



dr 



2r„ 



with 



r*L e -*$r- r^e-onn( 7 m 2 r) 

•/2t T ./2to T 



2(e B ) 3 4(e B ) 3 1 

Q12 (s, 77) = , Q22 (s, 77) = 



Q32(S,?7) 



3S ' " 3 S 1 _ jy2 ' 

(e B ) 3 4(e B ) 3 1 



2tt 



In ( 7 m 2 T ) £« + ^m 2 ^ - ^ (2) 



3s 



3s 1 - ?7 2 



Note that the function Qi(s, 77) = Qio(s,77) + <2i2(s,r7) is 
the expansion of Qi (s, 77) up to quadratic terms in the 
external field. Obviously, for tq — > the first term in 
Eq. (jT3"]l is logarithmically divergent (~ lnr _1 ). This 
divergence will be "reabsorbed" in the course of charge 
renormalization. 



with £(°) = -1/2B 2 . Here 
5m 



2 _ 3a mg 



2tt 



In 



2"" dm 2 



1 



(78) 



71x15x0 / 6 



(79) 



is the correction to the square of the "bare" electron 

mass (m 2 = m 2 , + 5m 2 ) [6, 30 32], whereas the renormal- 

(2) 

ized two loop term £^ R corresponding to each eigenmode 



11 



reads 



,(2) 
J 1R 



5a mg d£$ a f°° dr^_ m 2 r 
12tt dm 2 16tt 3 



/ ^e^ln ( 7 m 2 r) 
Jo T 



s coth(s) 



sinh (s) 



- 2 



32tt 3 J 



dTe- m « 7 



x / d77{Qi(s,r ? )-Q 1 (s,r ? )-^i(s,7 ? )} ) (80) 
Jo 

(2) _ 5aom|gg^ _^o_ A 00 dr ?T f , . 
^ 2R " 24n 9m 2 , 32^3 J r 3 in tf™^ 



3s coth(s) + 



-4- 



2s 2 



Of) 



sinh (s) 
/•l 

x e" m ° r / dT? {Q 2 (s, 77) - q 2 (s, 77) - jF. 2 (s, 77)} , (81) 
Jo 



T ( 2 ) _ _ 15a mg 9£j£ n,, 

'3R — 



24tt dm 2 , 32 
s coth(s 



3s 2 „ 2s 2 



sinh (s) 



32tt 3 



dr 



x e" m o r / d7 y {Q 3 (s,ry)-Q 3 (s,ry)-^3(s,r?)}.(82) 
./o 



Inserting Eqs. (|76H78p into Eq. (|60j) . the twodoop 
correction to the Euler-Heisenberg Lagrangian is given 
by 



£( 2 > = -^£(°) r^e-^ + W^V + 4 2 ) 
4tt 2 7 2tq r dm 2 

with 4 a) = ELi-c^- 

Considering Eq. (|54]l and Eq. ([83]) , we obtain 



(83) 



£ EH - -2 (0) ^3 ( Lo P ) + £ R 



an (l) 



9m 2 ' ' 



(84) 



where 



Note that eB = eoBo. So, the variable s is an invariant 
under the renormalization. Clearly, m must be replaced 
by m wherever it appears as well as ao — > a. Keeping 
this in mind, 



^EH — Ajji 



1(0) , (i) 



-R 



(2) 
iR 



(87) 



with 



a(2) _ 
-^lR — 



5am 4 b 



n 4 b f 00 ds s/b / 

-5- / "^e /b SCOth S 
7T 3 7 S^ V 



96 

,4u2 /-oo 



1 



16tt 3 7o s 3 V 7 bJ 



s coth(s) 



32tt 3 



sinh 2 (s) 

'a-'wfS^H- 1 - 



dse- s/b 



am b 



4u2 /•<» 



ds 



-s/b ln 



3s coth(s) + 



2s 



21 



sinh 2 (s) 



am b 



4-u2 r°° 



32tt 



3 



ds 



x e 



a(2) 
^3R 



/b / dr ? {Q 2 (s,r ? )-Q 2 (s ! r ? )-^ 2 (s ! r ? )}, (89) 



^e" s / b fscoth( B )-l-V 
s^ \ 3 



10 

15am b 

1927T 3 

,^4u2 r°° 



32 /i - jq 
s coth(s) 



3s 2 



sinh (s) 



-4- 



2s 2 
3 



32tt 3 



ds 



s/b / dtj{ 63(8,7?) -QaM-J^s, V)} (90) 



Z" 1 

^3{2loop) 



_ lim f°^V m ° T .(85) 



For T ^ the integral ^ e - m ° T = ln^ing-To) -1 
and therefore 



z- 1 

^3(2loop) 



1 + ^ln 



1 



^ln 



1 



t ->-o 37r ^mpTo/ 47r 2 \27mQro 

Now, we are able to identify 

d£ (1) 



4 1) (m 2 )^<(m 2 ) + f 5m 2 



dm 2 . 



(86) 



Sm 2 =0 



where m is the renormalized electron mass. The renor- 
malized charge and field strength are introduced by 
means of the relations e = eo-Z^^fioop) and B = 
B o-2 3( 2i / o 2 op)- Under this condition £^ 0) = £(°> Z7 1 



'3(2loop)' 



where we have replaced the integration variable r by 
s, and b = B/B c . Here we have used the function 

f»(s, 77) = |fii(s,r?) - Qi(s, 77) -^i(s,r?)| with Qi(s,r?) = 

Qi(s,77)/(eB) 3 , Qi(s,77) = Qi( S> 7?)/(eB) 3 and ^(3,77) = 
J0 i (s,r ? )/(eB) 3 . 



C. Asymptotic behavior at large magnetic field 
strength 

In the asymptotic region of superstrong magnetic field 
b > 1, Eq. flSBJ) behaves like (Ha 



*g>0» = (£)+>». , 



where 67r" 2 C'(2) = -0.5699610... with ((x) the Rie- 
mann zeta-function. 
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For the same magnetic field regime the leading term of 



( 2) 

£^ R is derived in Appendix [B] 

(2) _ am 4 b 2 



»(2) 

J 2R 



am 4 b 2 
32tt 3 



"1R 



AA In 



16?r 3 

b 

77T 



A/i, 



3 ln 



b 

77T 



.4 



(92) 



10 < b < 10 2 • 37r/a. Within this approximation, we 
find that the vacuum reacts paramagnetically and has a 
nonlinear dependence on the external field (see Eq. (f96|)). 

The two-loop correction is given by M^ 2 ) = Xa=i Mp" 1 

where Mj 2 '' = 9£;^/9B is the contribution correspond- 
ing to a photon propagation mode. Making use of Eqs. 
we find 



i(3) 
^3R 



32tt 3 



■A/a In — 



b 

") 7T 



3 ln 



b 

77T 



— \+B 



These expressions are calculated with accuracy of terms 
decreasing with b like ~ b _1 m(b) and faster. Here the 

numerical constants are A/i = 1.25, A/2 = h — 4< l2 , A/3 



_L 

3 ~~ P~ 

I + ^r 1 , A = 4.21 and B = 0.69. Note Af 2 + Af 3 = 1. 

Taking all this into account, the asymptotic behavior 
of the full two-loop term is 



x(2) _ 



£4? 



32tt 3 



77T 



2.4 



(93) 



which coincides with the results reported in references 

MM 



VI. VACUUM MAGNETIC PROPERTIES IN A 
SUPERSTRONG MAGNETIC FIELD 

A. Role of the photon polarization modes on the 
vacuum magnetization 

In presence of an external magnetic field, the zero- 
point vacuum energy <? vac is modified by the interaction 
between B and the virtual QED-particles. The latter 
is determined by the effective potential coming from the 
quantum-corrections to the Maxwell Lagrangian which is 
also contained within the finite temperature formalism. 
According to Eq. (|87|) it is expressed as 



<5\-. 



(2) 
iR 



(94) 



Consequently the vacuum acquires a non trivial magneti- 
zation M vac = — Sd'vac/c'B induced by the external mag- 
netic field. In what follows we will write 



M v 



= M« 

vac 



M ( 



2) 
vac 



(95) 



in correspondence with the loop-term £p . In this sense, 
the one-loop contribution at very large magnetic field 
b ^> 1 can be computed by means of Eq. (|9"Tj) and gives: 
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24tt 2 B c 
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(96) 



The dependence on the external field is shown in FIG. 
m The data depicted is obtained in the field interval 
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(99) 



According to these results, in a superstrong magnetic 

(2) (2) 
field limit, M^ ; < and M, ' < behave diamagnet- 

(2) 

ically whereas M3 > is purely paramagnetic (see 
Fig. [5]). Within the range of magnetic field values for 
which the photon anomalous magnetic moment is defined 
(10 14 G < |B| < 10 15 G) the vacuum magnetization den- 
sity is M 2 2) 10 9 erg/(cm 3 G). Moreover, while M^ 2) 

depends linearly on b, the contributions of the second 
and third propagation mode depend logarithmically on 
the external field. Note that the leading behavior of the 
complete two-loop contribution is 



M< 2 ) 



am b 
327r 3 B r 
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(100) 



which points out a dominance of the third mode. 
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FIG. 4: One-loop contribution to the vacuum magnetization 
density with regard to the external field. 



As it was expected M^/M^ 2 ' ~ a 1 . This ratio is also 
manifested between the corresponding magnetic suscep- 
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FIG. 5: Contribution of the vacuum polarization eigenvalues 
to the vacuum magnetization density with regard to the exter- 
nal field. Here 10 < b < 10 2 Sir/a. The dashed line represents 
the complete two-loop contribution to the vacuum. 
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FIG. 6: Contribution of the vacuum polarization eigenvalues 
to the vacuum transverse pressure density with regard to the 
external field. Here 10 < b < 10 2 • 3iv/a. The dashed line 
represents the complete two-loop contribution. 



tibilities (X« = dM^/dB). Note that 
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For magnetic fields b ~ 10 5 corresponding to |B| ~ 
10 18 G, the magnetic susceptibility reaches values the or- 
der of xW ~ 10 _4 erg/(cm 3 G 2 ) which exceeds the val- 
ues of many laboratory materials, for example Aluminum 
(X A i = 2.2 ■ lCr 5 erg/(cm 3 G 2 )). 



B. Transverse pressures 

Because of the anisotropy generated by B a magne- 
tized vacuum exerts two different pressure components 
[13, Hl| . One of them is positive (P|| = — <o vac ) and along 
B, whereas the remaining is transverse to the external 
field direction (Pj_ = -<f vac - M|B|). For b ~ 1 the 
latter acquires negative values. 

At very large magnetic fields (b > 1) the one- loop 
approximation of Pj_ can be computed by making use of 
Eq. flST} and Eq. ®. In fact 
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Therefore, at asymptotically large values of the external 
field, the interaction between B and the virtual electron- 
positron pairs generates a negative pressure which would 
tend to shrink inserted matter in the plane transverse to 
B. 

Again, the two-loop contribution can be written as the 
sum of the corresponding terms due to the vacuum po- 
larization modes 
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According to Eqs. (|92j) and Eqs. (|99|) they read: 
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with the complete two-loop term given by 
am 4 b 2 
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(107) 



For b ~ 10 5 , corresponding to magnetic fields B ~ 10 18 G, 
the transverse pressure generated by the first and sec- 
ond polarization mode is positive and reaches values of 
the order ~ 10 30 dyn/cm 2 and ~ 10 31 dyn/cm 2 , respec- 
tively (see Fig. [5]). In contrast, the contribution given by 
the third mode is negative with P^j ~ — 10 31 dyn/cm 2 . 
In the same context Pj 2 ' ~ — 10 31 dyn/cm 2 . The one- 
loop contribution is correspondingly even of the order of 



>(i) 



10 33 dyn/cm 2 



VII. CONCLUSION 

In the first part of this work we showed that an in- 
frared photon propagating in a strongly magnetized vac- 
uum (10 < b <C 37r/a), seems to exhibit a nonzero vector 
anomalous magnetic moment [see Eq. (|42|)]. We have 
pointed out that this quantity is a signal of the Lorentz 
symmetry breaking due to the presence of an external 
magnetic field. In addition, we have shown that m 7 
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arises due to the interaction between virtual electron- 
positron quantum pairs with the external magnetic field 
and its existence is closely related to the gauge invari- 
ance. In this context, we showed that m 7 can be decom- 
posed into two orthogonal components in correspondence 
with the cylindrical symmetry imposed by the external 
field. These components are opposite in sign, and only 
the one along B is conserved. We remarked that the 
photon paramagnetism (analyzed in [17]) is only associ- 
ated with a kind of photon magnetization rather than a 
magnetic moment. 

In the last sections of this work we discussed the ef- 
fect of the vacuum polarization tensor, which modifies 
the zero-point energy of the vacuum even in the absence 
of electromagnetic waves. In the limit of a superstrong 
magnetic field, the two-loop contribution of the magne- 
tization density corresponding to the second and third 
propagation mode depends nonlinearly on the external 
magnetic field and their behavior is diamagnetic and 
paramagnetic, respectively. On the other hand, the con- 
tribution coming from the first mode is diamagnetic and 
depends linear on B. We have seen that for very large 
magnetic field the contribution due to the third mode 
strongly dominate the analyzed quantities. In this mag- 
netic field regime the latter exerts a negative transverse 
pressures to the external field. On the contrary those 
contributions coming from the first and second virtual 
mode are positive. 

We want to point out that, although the decomposi- 
tion of the two-loop term in the Euler-Heisenberg La- 
grangian was considered in a magnetic background it re- 
mains valid also for the electric case. This fact should 
allow to study the role of the vacuum polarization modes 
in electron-positron production in a superstrong electric 
field. A detailed analysis of this issue will be presented 
in a forthcoming work. 



and its substitution into Eq. (1351) leads to: 
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Defining the photon magnetization by means of the rela- 
tion /i, 7 = — dui/dB we find 
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Note that for magnetic field strength 10 < b <C 37r/a 
one should treat g <C 1. As consequence Eq. (|A4|) be- 
comes 

iti = — g— ni cos</> and u\ = q— — nil sin^. (A5) 
' mo 2mo 



1. Relation between /i 7 and m 1 

According to these expressions we can write 



dB 



B. 



(A6) 



The second term in Eq. (|A6|) is a vector orthogonal to 
the external field. It can be expressed as 
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Appendix A: 



The term Z2 can be written as 
1 



Z2 



2A 



[B x k]< 



(B-k) s 
2A 



(Al) 



"liS 1 B = 0^— n_LCOs^ = — (n-ny) (n y x s 7 ) (A7) 
oB 2m 2m 

By substituting Eq. (g2j and Eq. ([AT]) in Eq. (|A6j) we 
obtain 



^ =0 2nT [ n + ( n ' n ll) n ll] X 



7 v o , , . L V 11/ II J 7" 

(A8) 

Substituting Eq. (KB in Eq. (gTJ gives 



m 1 B = —fi B 



2m 



cos <pn ± ■ B (A9) 



Because of the fact that the second term on the right- 
hand side is projected out by the scalar product (nj_ B = 
0), the physical consequences of fi^ might be analyzed by 
means of m 7 . So, in a good approximation 



(A10) 



the remaining transversal component involved in m 7 is 
not relevant because is projected out by the scalar prod- 



uct m 7 • B. 
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2. fA for b > 3ir/a 



Appendix B: 



Note, in addition, that for purely perpendicular prop- 
agation, Eq. (|A4I) reduces to: 



/4=fl^/[kJ(l + 0(b)) 
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In particular for b — > oo 
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The aim of this appendix is to find the asymptotic be- 

(2) 

havior of £V R y for b>l. For this purpose, it is convenient 
to write 
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The second term in Eq. (|B1[) is defined as: 

4i 2 rl p oo 

£ P = -^^ G *( b ) with G i( b )= / df ? / dsc-^Ms,^) 
■i^ Jo Jo 

with fi(s, 57) = |Qi(s,r?) - Qi(s, 17) - ^i(s,7?)| . 
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1. Leading behavior of £J R in an asymptotically large magnetic field 

In order to determine the leading asymptotic- magnetic field term of £^ , we substitute Eq. (|91j) into Eqs. (|B2|) . 
which gives: 
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(B5) 



Note that £2 and £3 can be expressed as 



with 



Si + 2S, £ 3 = £i-2E-2b 
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To derive the second expression in Eq. (|B6[) we have used 
the identity 
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and an integration by parts. 

Note that Si converges even without the exponential 
factor which approaches to 1 for b — > oo. By using 
MATHEMATICA code we find Si ~ 0.19. 

S does not involve singularities in the integrands at 
s = 0, but would diverge at s — > oo if one sets the limit- 
ing value exp(— s/b) = 1. For that reason we divide the 
integration domain into two parts: 
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and T an arbitrary positive number. 

Now, we can omit the exponential in Sl since the re- 
sulting integral converges anyway: 
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Substantial simplification is achieved by splitting the 
integrand of Sh into its parts and neglecting the expo- 
nential factor exp(— s/b) whenever is possible 
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For s — > oo the leading term of coth(s) ~ 1. Having this 
in mind, we compute the integrals: 
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In the latter we have neglected terms decreasing as ~ 
T/b. Numerical calculation using MATHEMATICA code 
gives T ~ 7r by minimizing S (dS/dT = 0). Therefore, 
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and according to Eq. (|B6 
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The substitution of the latter into Eqs. (|B5|) yields 
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2. Leading behavior of £^ R in an asymptotically 
large magnetic field 

(2)2 

The asymptotic behavior of £,- R ; is obtained from the 
integral Gi. We start our analysis by dividing the inte- 
gration domain into two regions: 

/•T /*1 roo /*1 

Gi= ds / dr]...+ ds / d?7 . . . (B17) 

Jo JO JT JO 

with T > 0. We denote the corresponding integrals by 
Gj and G^ , respectively. In the latter we replace the 
upper integration limit over 77— variable by a parameter 
770 = 1 — T/s > 0. In order to find the behavior of G?^ 
we write 



G[ L) = ijf ds^ d^f i (s,ry)exp(-s/b) (B18) 
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where we have set exp(s/b) ~ 1 since fi converges within 
the corresponding integration domain. For T — > the 
behavior of the remaining integrand is 
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The integrals over s and 77 are trivial to perform and give: 
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Let us consider now, the contributions coming from 
the second integration domain which concerns to G- . 
In order to this we consider the asymptotic expression 
of fi(s, 77) for s — > 00. First of all, the asymptotic ex- 
pansion of Eqs. ([5]) in powers of exp(— s) and exp(s?7) 
produces an expansion of Eq. in a sum of contribu- 
tions coming from the thresholds, the singular behavior 
in the threshold points originating from the divergencies 
of the s— integration in Eq. Q near s = 00 as it was 
developed in [9j. The leading terms in the expansion of 
Eqs. iJSJ) at s — > 00 are 
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These expressions are in correspondence with the lowest 
threshold (n = 0, nf = 1 or viceversa for i = l,n = n' = l 
for i = 3, and n = n' = for % = 2). The fact that 
M(oo, 77) ~ 1/2 and N(oo, rj) 3> M(oo, rj) allow us to write 
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with V(s, 77) + W(s, ry) ~ 4N _1 (s, rj). Considering only the 
terms which decrease most slowly as a function of s, we 
find 
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Additionally, the most significant terms arising from Q; 
S)i for s — > 00 are 
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The behavior of fi — exp(— b/s)(Qi — Qj — i5i) for s — > 00 



is find out by considering Eqs. (|B24HB25l) . In fact 
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Except for the last term in fi^^ and fal^^ , we may 
use the exp(— s/b) ~ 1 so that 
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Performing the integration over 77, G^ are given by 
f°° p~ 2s 

G 1 H) ~ 4 / ds [Ei(2s - T) - Ei(T)] (B28) 
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where Ei(— T/b) is the exponential-integral function 
whose asymptotic expansion for very large magnetic field 
b — > 00 is 



ds / s 
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This expression is calculated with accuracy of terms that 
decrease with b. To be consistent the last integral of 
Eq. (|B28|) and Eq. (|B30[) can be neglected as well, since 
the terms decrease as fast as ~ b _1 lnb and ~ b _1 . The 
remaining integrals present in Eq. (|B29|) and Eq. (IB30[) 
depend on the parameter T. Taking all this into account, 
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the leading asymptotic behavior of Gi for b — ¥ oo reads and 
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where we have used Eq. (|B31|) . Here, the numerical 
constants Ci,3 are determined by imposing the condition 
dCi, 3 /dT = 6 with 



C = *J T ^ [Ei(2s-T)-Ei(T)]-^ 



ds In 



/2s \ e" 2s 



















,(B33) 









:?-o 







41T 2 
1620' 



This yields 



Ci w 2.67, and C 3 « 0.18. 
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Note that there is not value fulfilling the condition 
dC2/dT = 0. However, in order to compute it, we first 
set ^ £i=i Gi = 0. This condition leads to T ~ 0.46 
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Obviously, G2 = G — Gi — G3. Thus, by taking into 
account Eqs. (|B32j) and Eq. (|B35[) we obtain C 2 « -4.68. 
Such that 
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Substitution of Eqs. (|B37[) in Eq. (fB4l) allows to obtain 
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C 2 with C 2 = -4.68, (B39) 
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Finally, the behavior of £^ presented in Eq. is 
derived by inserting the expression bellow Eqs. (|B16|) 
and Eqs. (|B38HB40j) into Eq. (|BT]l . 
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